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ABSTRACT 

The G305 H ii complex (G305. 4+0.1) is one of the most massive star forming structures 
yet identified within the Galaxy. It is host to many massive stars at all stages of 
formation and evolution, from embedded molecular cores to post main-sequence stars. 
Here, we present a detailed near-infrared analysis of the two central star clusters 
Danks 1 and Danks 2, using HST-hNICMOS imaging and VLT-MSAAC spectroscopy. 
We find that the spectro-photometric distance to the clusters is consistent with the 
kinematic distance to the G305 complex, an average of all measurements giving a 
distance of 3.8±0.6kpc. From analysis of the stellar populations and the pre-main- 
sequence stars we find that Danks 2 is the elder of the two clusters, with an age of 
S^iMyr. Danks 1 is clearly younger with an age of I.SIq jMyr, and is dominated by 
three very luminous H-rich Wolf-Rayet stars which may have masses ^ lOOM©. The 
two clusters have mass functions consistent with the Salpeter slope, and total cluster 
masses of 8OOO±15OOM0 and 3000±800Mq for Danks 1 and Danks 2 respectively 
Danks 1 is significantly the more compact cluster of the two, and is one of the densest 
clusters in the Galaxy with log(/9/M0pc~'^) = 5.51q^4. In addition to the clusters, 
there is a population of apparently isolated Wolf-Rayet stars within the molecular 
cloud's cavity. Our results suggest that the star-forming history of G305 began with 
the formation of Danks 2, and subsequently Danks 1, with the origin of the diffuse 
evolved population currently uncertain. Together, the massive stars at the centre of 
the G305 region appear to be clearing away what is left of the natal cloud, triggering 
a further generation of star formation at the cloud's periphery. 

Key words: (Galaxy:) open clusters and associations: general (Galaxy:) open clus- 
ters and associations: individuahDanks 1, Danks 2 stars: Wolf-Rayet stars: formation 
(ISM:) H ii regions ISM: clouds 



1 INTRODUCTION 

Massive stars have a profound effect on their wider Galac- 
tic environment, via the production of copious quantities 
of ionising radiation, and from the input of mechancial en- 
ergy and chemically processed matter into the interstellar 
medium (ISM). For these reasons, an understanding of their 



lifecycle is of importance to many areas of astronomy. Unfor- 
tunately, a number of questions regarding this still remain 
unanswered, with the nature of their formation mechan- 
sism(s) being particularly opaque. While growing observa- 
tional evidence suggests that stars between 20-40Mq may 
form via disc mediated accretion - in a manner analagous 
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to their lower mass counterparts (e.g. W33A & W51N, see 
iDavies et al.l I2OI0I : IZapata et al. 20091) - it is still not clear 
how more massive stars form teinne cker fc YorkdlioOTl ): de- 
spite compelling observational evidence for stars with masses 
signifi cantly in excess of AQMp, (WR20a, NGC 3603- Al, and 
R145. iBonanos et al ] |2004l : lRauw et al.ll2005l : ISchnurr et alJ 
I2OO8I . I2OO9I ). 

Stellar hierarchies appear to be a signature of star for- 
mation, with stars predominantly forming in cluste rs and, 
in tu rn, clusters forming in larger complexes (e.g. iLarsenl 
I2OO4I ). Such structures, spanning tens to hundreds of par- 
sees, are most readily identifiable in exter nal star-forming 
galaxies such as M51 jBastian et all 120051 ). While ages for 
individual clusters within the M51 complexes are difficult to 
determine, it appears that they are likely to be rather youth- 
ful (e.g. <10 Myr) and massive (S-SOxlO^M©). Similar com- 
plexes are seen within interacting, starbursting galaxies such 
as the Antennae, where star formation rates are an order of 
a magnitude higher than in M5 1, in turn yielding individ ual 
clusters with masses >1O'^M0 (jBastian et al.ll2005l . [200^ ). 

A precise understanding of the nature of such complexes 
would be invaluable for the following reasons: (i) they appear 
to represent a ubiquitous mode of star formation in star- 
burst galaxies; and (ii) by virtue of their masses they pro- 
vide a statistically well-sampled stellar mass function. Un- 
fortunately, the distances of their host galaxies and compact 
nature conspire to make the determination of the properties 
of individual clusters - let alone stars - observationally chal- 
lenging. Therefore, one might ask whether such structures 
are present within our own Galaxy? The recent detection 
of a number of massive (>10''Mq) Red Supergiant (RSG) 
dominated clusters at the b ase of the Sc utum-Cru x arm is 
suggestive of such a complex (iFiger et al.1[20 06: Da vies et al.l 
I2OO7I : IClark et al.ll2009bl : iNegueruela et ahfeoig . I2011al ). al- 
though their spatial extent (~100pc) and age spread (~10- 
20 Myr) currently preclude an unambiguous association with 
a single, physically distinct structure. 

One observational approach to overcome such uncer- 
tainties is to search for young, massive clusters still embed- 
ded in their natal giant molecular cloud (GMC) and/or asso- 
ciated giant H 11 region. Such a strategy guarantees the youth 
of such a complex, potentially enabling individual examples 
of massive young stellar objects (MYSOs) to be identified, 
and ultimately the global star formation history from a spa- 
tially resolved census of the (proto-) stellar populations. The 
latter goal is particularly important, since the limited spa- 
tial resolution of such objects in external galaxies precludes 
a detailed analysis of the processes by which the GMC is 
converted into stars and star clusters. 

A number of GMCs which appear to contain both mas- 
sive (>10^ Mq), young clusters as well as deeply embed- 
ded MYSOs have been identified (e.g. W49A, W51 and the 
Carina nebula; IaIvcs & liomcicr 20031: iKumar et al.ir2004 



Table 1. Read-sequences and total integration times employed 
for each filter during the NICMOS observations. 



ISmith fc Brooks |2007|). Another such re gion is the G305 star 



forming complex (jClark fc Porterll2004l . ^=305.4, h=+G.l; ). 
Located in the Scutum-Crux arm at an estimated distance 
of ~ 4 kpc, it has the form of a large tri-lobed cavity with 
a maximum extent of ~34pc, delineated by both mid- and 
far-IR, sub-millimetre and radio emission and centred on the 
young clusters Banks 1 & 2 (Fig. [T|. A large contingent of 
massive stars is inferred from the ionising radiation required 
to support the total radio flux (equivalent to the output from 
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SAMP-SEQ 


NSAMP 


Tint (s) 


F160W 


STEP2 


17 


16S 


F222M 


STEPS 


12 


44S 


F187N 


STEPS 


10 


240 


F190N 


STEPS 


10 


240 



>30 canonical 07 V stars). In addition, there are numerous 
signposts of ongoing star formation in the cloud's periph- 
ery, in the form of deeply embedded MYS Os, compact Hll- 
regions, and meth anol and wa t er masers (UJrauhart et al.l 
2007, 2009; Hindso n et aLll2010l : IClark et al.ll201ll ). Finally. 
far-IR and sub-millimetre continuum observations reveal the 
presence of a significant reservoir of cold molecular mate- 
rial (>10^M(:)) availa ble to fuel further star-forming activity 
iHindson et alj|2010l . Clark et al. in prep). 

The overal morphology of the complex is strongly in- 
dicative of a number of epochs of sequential star forma- 
tion, initiated and sustained by the action of the two cen- 
tral clusters. In this paper, we present a near-IR analysis of 
these clusters, using high-resolution photometry and spec- 
troscopy, in order to determine their masses and ages, and 
consequently whether their properties are consistent with 
such an hypothesis. We discuss and present these data in 
Section 2, report our analysis in Section 3 and discuss these 
result both in the context of the G305 complex and also in 
comparison to other star forming regions in Section 4, before 
summarising our conclusions in Sect. 5. 



2 OBSERVATIONS & DATA REDUCTION 
2.1 Photometry 

2.1.1 Observations 

Images of the two clusters were obtained with 
HST/NICMOS on 16 July 2008. We used the NIC3 
camera, which has a field-of-view of 51.2" x 51.2" and a 
pixel scale of 0.2". The clusters were imaged through each 
of the filters F160W and F222M, as well as the narrow-band 
filters F187N and F190N which are centred on Pq and the 
neighbouring continuum respectively. In addition to the 
clusters, in order to characterize the foreground population 
we also imaged nearby control fields through the F160W 
and F222M filters. The observed fields are indicated in 

Fig.m 

In our observations we employed a spiral dither pat- 
tern with six separate pointings, each offset from the last 
by 5.07". By employing sub-pixel dithering we minimise the 
effects of non-uniform intra-pixel sensitivity. We used the 
MULTIACCUM read mode, using read-sequences and pat- 
terns that provided good sampling coverage over a large dy- 
namic range. The sampling sequences and total integration 
times we used for each filter are listed in Table [T] 

2.1.2 Data reduction and analysis 

Following the guidelines of the NICMOS Data Handbook 
v7.0 we performed the standard reduction steps of bias sub- 




Figure 1. Wide-field Spitzer/GLIMPSE image of the G305 complex at S.S^tm, with the NICMOS fields overplotted. The 5.8/im band 
traces the PAH emission in the complex, and therefore the interface between ionised and molecular gas. 



traction, dark-current correction and flat-fielding using the 
custom-written NICMOS software CALNICA. as our dithered 
observations made three sub-samples of each pixel we re- 
sampled each dithered observation onto a finer grid contain- 
ing a factor of 3 more pseudo-pixels. The six images were 
then mosaicked together, again using the custom-written 
NICMOS software. The final mosaics of each cluster are 
shown in Fig. [2] 

In extracting the photometry from the images, we ex- 
perimented with several methods and algorithms. We found 
that the best results were obtained (fewest spurious source 
detections, most effective doblonding of sources) using the 
STARFINDER routines (Diolaiti et al. 2000) which run under 
the IDL environment. The algorithm uses the image point- 
spread function (PSF) in order to locate point-sources. We 
found that the algorithm was most effective when using the 
synthetic PSF generated by tinytim. To check for consis- 
tency between this method and regular aperture photome- 
try, we ran both STARFINDER and daophot on the relatively 
uncrowded control fields and found no systematic differences 
between the two. 

In order to characterize the statistical completeness in 
each field observed, as well as make empirical measurements 
of the photometric errors, we performed artificial star exper- 
iments on both cluster fields. We used tinytim to generate 



the artificial stars with a luminosity function identical to 
that observed in each field. We then randomly sampled arti- 
ficial stars from this luminosity function and added them to 
the image. No more than 100 stars were added to the image 
so as not to significantly alter the level of crowding. The pho- 
tometry algorithm was then re-run on the artificial image, 
and the output photometry and astrometry compared with 
the locations and magnitudes of the input stars. For a star to 
be considered 'recovered' we specified a maximum distance 
separation between input position and detected position of 
0.22", or roughly la of the PSF. Additionally, if the out- 
put magnitude of a star was more than a factor of 2 greater 
than the input magnitude then it was assumed that that 
star had been blended with a brighter one, and that the in- 
put star was lost. The 50% completeness limits were found 
to be at mieow = 18.8 and m222M = 17.4 for Danks 1, and 
T7^l60w = 19.2 and m222M = 18.3 for Danks 2. The slightly 
fainter detection limit for Danks 2 with respect to Danks 1 
is due to the reduced level of crowding in that field. 
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Figure 2. F160W mosaics of eacli cluster; left: Danks 1, right: Danks 2. In each image nortli is up and east is left. The stellar identifications 
(see Table [2]| are denoted by the red circles. 



2.2 Spectroscopy 

2.2.1 Observations 

Spectra were obtained of several stars in each cluster using 
the ESO ISAAC instrument mounted on UTl of the VLlQ. 
Observations were taken in good weather on the nights of 
June 27th to July 2nd 2006. The instrument was used in 
medium resolution mode with the 0.3" slit, and two over- 
lapping wavelength settings were used per target to achieve 
a spectral range of 2.04-2.21/im at a spectral resolution of 
A/AA=8900. 

The observing strategy was to align the slit in order to 
obtain spectra of two programme stars simultaneously. The 
stars were nodded along the slit in an ABBA pattern. Inte- 
gration times per slit position were limited to '^lOO seconds 
to allow for accurate subtraction of the sky emission, with 
each integration split into 4-10 seperate read-outs to avoid 
saturation around strong stellar emission lines. The total in- 
tegration times per star were between 400-500 seconds. The 
star HD 113457 (spectral type AOv) was observed after ev- 
ery two target observations as a measure of the atmospheric 
absorption. Continuum lamp exposures were taken just as 
frequently in order to correct for any fringing on the detec- 
tor. 

2.2.2 Data reduction 

Nod pairs were subtracted from one another to remove bias 
level, dark current and sky emission lines. Before extracting 
the spectra, the frames were corrected for the degree of warp- 
ing which is present in ISAAC spectroscopic data. Warping 
in the spatial direction was characterized by fitting polyno- 
mials to the stellar spectral traces across the detector. To 

1 ESO programme ID 077.C-0207(B), PI J.S. Clark 



measure the warping in the dispersion direction a 'sky lines' 
image was created by summing all nod pairs of the same tar- 
get and subtracting the stars. The sky lines were then fitted 
with polynomials to wavelength-calibrate the data and to 
measure the degree of warping in the dispersion direction. 
Using the fits to the stellar traces and the sky lines each 
frame was then resampled onto a linear, wavelength cali- 
brated grid. The wavelength solution of each frame had an 
r.m.s. of 0.2-0.3 pixels, or 3-6 km s~^. 

In some observations, it was noticed that a degree of 
diffuse line-emission was present due to ionized nebulae in 
the vicinity of the clusters. This emission was subtracted 
from each frame by fitting gaussian profiles to the spatial 
variations and interpolating across any stellar traces that it 
intersected. 

Spectra were extracted by summing across the pixels 
around each stellar trace. Cosmic ray hits and bad pixels 
were rejected by comparing repeat observations of the same 
star. Before dividing through by the telluric standard spec- 
trum, we first fitted the standard's Br 7 absorption feature 
using a Voigt profile, and corrected for the continuum slope 
by dividing through by a black-body appropriate for the 
star's spectral type. The standard star spectrum was then 
cross-correlated with the science target to correct for any 
small sub-pixel shifts which may produce artifacts in the fi- 
nal spectrum. The signal-to-noise ratio (SNR) of the fully 
reduced target spectra was estimated from flat regions of 
continuum, and was found to be typically 150-400. The fi- 
nal spectra for the stars in Danks 1 and 2 are shown in Fig. 
[3] and Fig. [4] respectively. 

2.2.3 Spectral classification 

With the exception of one object (D2-2), we found that all 
stars observed exhibited spectral features associated with 
hot stars. The diagnostic features observed were Br 7, He I 
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Figure 3. Spectra o f the stars in Dank s 1. Also shown arc compariso n spectra of template stars taken from I Hanson et al.l ([199^, [ioO^ , 
iMartins et al.l ll2008l ). rFieer et al.l lll997l ). and lWallace fc Hinkld l l 19971 ). The locations of key diagnostic spectral lines have been indicated. 
The three panels show the WNh stars {top), the early-mid O stars (middle), and the late-O / early-B stars [bottom). 
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Figure 4. Same as Fig. [3] but for Danks 2. Tlie tliree panels show the OB stars (top), the WCL star (middle), and the yellow star 
{bottom). 
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(2.056/im, 2.112^im, 2.165/im), Hen (2.189^im), the blend of 
Cm, Niii and Oiii at 2.115^m, Niii (2.103/im), and Civ 
(2.070Mm, 2.079Mm). 

In orde r to classify the spectra o f the stars, we use d 
the works o f lHanson et all (Il996l lioOsh.lMorris et alj j 19961 ). 
iFiger et all l|l997h and lBohannan fc Crowtheij l| 19991 ). Also, 
since a number of Danks stars have ii'-band spectra which 
are remarkably similar to those in the Arches cluster, we 
m ake comparisons to t he spectra of those stars presented 
in iMartins et al.l l|2008h . Given the comparative paucity of 
emission features suitable for classification, spectral types 
determined solely from K band spectra are inevitably less 
precise than those obtained from optical data. Nevertheless, 
our spectral classification scheme is as follows; 

• WNLh: broad emission lines of Br7, the 2.115/im com- 
plex. Civ and Niii. The Hei (2.056^im) and Heii(2.189^im) 
features have complex structures and/or P Cygni-type pro- 
files. Emission from C iv is either absent or very weak. 

• O hypergiant (Olf): broad emission of Br 7 but with 
low contrast to the continuum; 2.115pim complex and Civ 
in emission; He I (2.056/im, 2.112/im) and Hen (2.189/im) 
in absorption. 

• 0-type: broad Br 7 in absorption, though for super- 
giants it may be in emission. C iv emission and He 11 
(2.189/im) absorption is seen in spectral types no later than 
07 and 08 respectively. Hel (2.112/im) is in absorption with 
the Hel 2.113/im multiplet in emission for types later than 
08. Hel (2.056/im) is heavily dependent on the wind prop- 
erties and may be in either emission or absorption. 

• late 0-/early B-type: Br7, Hel (2.112/im) and Hel 
(2.162/im) in absorption. Absence of Civ in spectral-types 
later than 08, and no Hell types later than 09. The Hel 
(2.112/im) feature is absent in spectral-types later than 
~B5, while Hel (2.056/im) is not seen beyond ~B3 - though 
we note again that this line is very sensitive to wind density. 

The other spectra which do not fall into these categories 
are easily classifiable: D2-3 has the broad CiV and 2.115/im 
emission lines of a WC7/WC8 star (see Fig. [J] centre panel), 
and has been previously identified bv Mauerh an et al.. (,200 9, 
classified by those authors as WC8); while D2-2 has the 
spectral features of a late F-type star but without the dense 
molecular absorption of an early G-type star (Fig.|4l bottom 
panel). The near-IR colour of this star suggests that it may 
be a foreground star. 

The spectral types attributed to each star, along with 
their coordinates and photometry, are listed in Table [2] In 
Danks 1, we found three stars with spectral features at- 
tributable to H-rich WR stars (type WNLh), stars Dl-1, 
Dl-2 and Dl-fl Stars Dl-4 and Dl-6 have strong Civ in 
emission as well as very broad Br 7 emission, while Hel 
(2.112/im) is in absorption. These stars are assigned spec- 
tral types O6-81B Stars Dl-7, Dl-8, Dl-9 have weak Civ 
and Hell, with no Hel (2.112/^m), and are classified as 04-6. 



2 iMauerhan et al.l 1 I2OIII) also detected Dl-1 and Dl-5. They clas- 
sified Dl-5 as WN9h, consistent with our classification of WNLh. 
However, Dl-5 was classified as WN9, whereas our high signal- 
to-noise spectrum detected weak emission from Br 7 and so was 
classified WNLh^ 

3 We note that iMartins et al.l ||2008| '1 assigned types of 04-6If-|- 
for stars with similar spectral appearance in the iC-band. This is 



Dl-10 is similar to these stars except that Hel (2.112/^m) 
is in absorption, implying a slightly later spectral type of 
06-8. Stars Dl-11 and Dl-12 have no Hell and weak Hel 
2.112/im, and are classified as 08-B3. The width of the ab- 
sorption lines in the spectra of Dl-10, Dl-11 and Dl-12 
mean that they are likely to have luminosity classes of v- 
III. Star Dl-3 has very narrow absorption lines of Br 7 and 
Hel (2.112,2.162/^m) while Hel 2.056/^m is in emission, and 
so is considered to be a 08-B3 supergiant. 

In Danks 2, aside from the WC and yellow star, we find 
3 stars with very weak C IV emission, weak Br 7 absorption 
and Hel (2.112^m) absorption (D2-5, D2-8, D2-9), which 
we classify as 06-8. A further three stars (D2-4, D2-6, D2- 
7) are spectroscopically similar, but have no Civ, and so 
have slightly later types of 08-9. The remaining stars (D2- 
1 and D2-7) are classified as 08-B3, with the former star 
having the narrower lines and Hel 2.056/im emission typical 
of a supergiant. 

Finally, we remark that from analysis of the P a fiuxes 
of the stars in each field (see Sect.fS]), it is unlikely that there 
are any further strong emission-line stars in either cluster, 
aside from those presented here. 



3 RESULTS & ANALYSIS 

Our data have revealed a large number of massive stars 
within the two central clusters. These are in addition to 
the nu merous and apparently isolated massive stars within 
G305 (|Shara et al.ll2009l : IMauerhan et al.ll201ll ). In Table [2] 
we list the astrometry, photometry and spectral types of 
these stars. For completeness, in Table [3] we list the other 
known massive stars in the vicinity of the two clusters. 

In Figs [S] and [5] we plot the results of the photometry 
for each cluster. In the left panel of each figure we show 
the colour-magnitude diagrams (CMDs) of the cluster and 
control fields. In the centre panel we have decontaminated 
the cluster field of foreground stars using the control field: 
for every star in the control field, we subtract a correspond- 
ing star in the cluster field with /222M and (/leow — /222M) 
within the photometric errors at that brightness. As the for- 
mal photometric errors can be very small for bright objects, 
we specify a minimum 'errorbox' of size /222M = iO.lmags, 
(/16OW — /222m) = ±0.14mags for the decontamination al- 
gorithm. In the centre panels we also illustrate the bright- 
nesses and colours of zero-age main sequence (ZAMS) stars 
at the distances of the clusters (see Sect. [372] for discussion 
on cluster distances In the right-hand panel we show the 
photometry across P a. 

In the following sections we use the information in these 
plots, as well as the spectral types of the stars, to determine 
the physical properties of each cluster. In our analysis we 



representitive of the uncertainty when classifying hot stars in this 
spectral window. 

* These ZAMS tracks were computed by first taking the 
masses, temperatures a nd luminosities of ZAMS stars from 
iMevnet fc MaedeJ | |200C| ) . T he relation between sp ectral type and 

Martins et al.l ||2005| ') for stars with 



temperature was taken from 



masses >15M0 and I Johnson for the rest. Infrared mag- 

nitudes and colou rs were taken from [Martins fc" Pled | |2006|'| and 
iKoornneej l ll983l) . 
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Figure 6. Same as Fig. [5] but for the field centred on 



Danks 2. 
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Table 2. Coordinates (J2000), spectral types and photometry of the stars with spectroscopic data. ^This photometry is taken from the 
2MASS catalogue H-hanA data, as our F160W data was corrupted for this star. 



Star 


Alt. ID 




RA 






DEC 




Spec Type 


"IF160 


mF222 


mFi87 


nipigo 


Danks 1 


























Dl-1 


MDM 8 


13'* 


12™ 


28.49= 


-62° 


41' 


43.46" 


WNLh 


7.260t 


6.620 


6.094 


7.002 


Dl-2 




13 


12 


24.95 


-62 


41 


59.92 


WNLh 


8.158 


7.460 


6.941 


7.752 


Dl-3 




13 


12 


23.69 


-62 


42 


0.99 


08-B3I 


8.223 


7.618 


7.840 


7.925 


Dl-4 




13 


12 


26.16 


-62 


41 


57.50 


06-8If 


8.934 


8.262 


8.166 


8.511 


Dl-5 


MDM 7 


13 


12 


28.47 


-62 


41 


50.72 


WNLh 


8.834 


8.306 


7.943 


8.518 


Dl-6 




13 


12 


26.22 


-62 


42 


9.37 


06-8If 


9.220 


8.647 


8.595 


8.857 


Dl-7 




13 


12 


26.80 


-62 


41 


56.36 


04-6 


10.161 


9.540 


9.685 


9.754 


Dl-8 




13 


12 


22.84 


-62 


41 


48.60 


04-6 


10.355 


9.685 


9.817 


9.905 


Dl-9 




13 


12 


26.02 


-62 


42 


15.59 


04-6 


10.806 


10.239 


10.368 


10.398 


Dl-10 




13 


12 


24.50 


-62 


42 


8.52 


06-8 


11.319 


10.612 


10.805 


10.864 


Dl-11 




13 


12 


26.32 


-62 


42 


5.78 


08-B3 


11.804 


11.255 


11.431 


11.455 


Dl-12 




13 


12 


25.62 


-62 


42 


5.13 


08-B3 


12.227 


11.598 


11.800 


11.836 


Danks 2 


























D2-1 




13 


12 


56.34 


-62 


40 


27.78 


08-B3I 


7.750 


6.810 


7.027 


7.124 


D2-2 




13 


12 


59.90 


-62 


40 


39.71 


F8-G1 


9.131 


8.912 


8.935 


8.997 


D2-3 




13 


12 


57.63 


-62 


40 


59.42 


WC7-8 


9.911 


9.146 


8.392 


9.443 


D2-4 




13 


12 


56.20 


-62 


40 


51.11 


08-9 


10.108 


9.537 


9.616 


9.783 


D2-5 




13 


12 


54.48 


-62 


41 


4.60 


06-8 


10.001 


9.605 


9.625 


9.717 


D2-6 




13 


12 


58.22 


-62 


40 


36.64 


08-9 


10.125 


9.631 


9.740 


9.813 


D2-7 




13 


12 


58.56 


-62 


40 


54.84 


08-9 


10.185 


9.669 


9.755 


9.871 


D2-8 




13 


12 


54.37 


-62 


40 


45.48 


06-8 


10.333 


9.830 


9.967 


10.034 


D2-9 




13 


12 


57.38 


-62 


40 


1.43 


06-8 











Table 3. Coordinates (J2000), spectr al types and 2MASS ph oto metry of the other c onfi rmed massive stars in the G305 complex, take n 
fro m the literature. Refere nces are: 1 - iMauerhan et al.l l l201lD : 2 - IShara et ahl | |2009|) : 3 - Ivan der Huchtl 1I2OOII '): 4 - iLeistra et al] l l2005f) : 
5 - IClark fc Ported |200l)- 



ID 


Alt. ID 




RA 






DEC 




Spec Type 


J 


H 


Ks 


Ref. 


MDM 3 




13'* 


12" 


09.05= 


-62° 


43' 


26.7" 


WN8-9 


10.21 


8.57 


7.58 


1 


S09 845-34 


MDM 4 


13 


12 


21.30 


-62 


40 


12.5 


WC8 


10.75 


9.57 


8.77 


1,2 


MDM 5 




13 


12 


25.46 


-62 


44 


41.7 


WN9 


9.81 


8.48 


7.65 


1 


S09 845-35 


MDM 6 


13 


12 


27.66 


-62 


44 


22.0 


WC7 


13.16 


11.82 


10.71 


1,2 


WR48a 




13 


12 


39.65 


-62 


42 


55.8 


WC6 


8.74 


6.80 


5.09 


3 


S09 847-8 




13 


12 


45.35 


-63 


05 


52.0 


WN6 


13.06 


11.34 


10.26 


2 


L05-A1 




13 


11 


41.04 


-62 


32 


50.8 


05-6 I 


11.75 


10.39 


9.58 


4 


L05-A2 




13 


11 


33.88 


-62 


33 


27.1 


BO-1 V 


12.31 


11.02 


10.34 


4 


L05-A3 




13 


11 


39.50 


-62 


33 


28.2 


B2-3 V 


14.06 


12.65 


11.97 


4 


MSX305.4013-f00.0170 




13 


13 


02.04 


-62 


45 


03.3 


WCL 


6.59 


5.00 


3.95 


5 



will equate the NICMOS filters F160W and F222M with 
the photometric bands H and K respectively. 



Danks 2 respectively 0. We adopt these values throughout 
the rest of this work. 



3.1 Extinction 

The extinction towards each cluster is derived from the 
brightest stars in the decontaminated CMDs. The median 
F160W-F222M colour is obtained from all stars in each clus- 
ter with F222M magnitudes brighter than 13 whose locus 
is close to the ZAMS, since these stars show no evidence 
for circumstellar extinction. As the intrinsic F160W-F222M 
colours of hot stars are all approximately zero, we take the 
median colour as a colour excess, and conver t it in to an 
extinction using the law of iRieke fc Lebofskvl (|l985l ). We 
find average line-of-sight extinctions for the two clusters of 
Ak = 1.11 ± 0.16 and Ak = 0.92 ± 0.29 for Danks 1 and 



3.2 Distances 

In order to calculate the distance to the clusters, we employ 
two independent and complementary methods. We use vari- 
ous radial velocities measurements for the G305 complex to 
determine a kinematic distance, and the photometry of stars 
in each cluster to calculate spectro-photometric distances. 



^ These extinctions were computed using a value of a = —1.53 for 
the slope of the IR extinction law. If a slope of a = —2.14 is used 
llStead fc H oarc 2009), we find Ak = 0.6±0.1 and Ak = 0.6±0.2 
for Danks 1 and Danks 2 respectively. 
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Figure 7. Galactic rotation curve in the direction of G305 (solid 
line). The dotted lines show the extremes obtained for the rotation 
curve when the uncertainties in Galactocentric distance and Solar 
rotation velocity are taken into account. 

3.2.1 Kinematic distance to G305 

Many Hll-regions and young stellar objects (YSOs) are 
found in the dense molecular material surrounding the two 
clusters, and for many of these objects radial v elocity mea- 
surem ents exist. A search of the 7?M5 database l|Hoare et al.l 
l2005l fl of YSOs in the region displayed i n Fig, [l] yielded 
15 o b jects with known radial velocities (jUrauhart et al.l 
12001 I2OO9I . and refs therein). The mean radial velocity 
and standard deviation is wlsr = —39.4 ± 3.0 km s~^. 
In Fig. [7] we compare this value to the Galactic rotation 
cu rve in the d irection of G305. We use the rotation curve 
of iBrand fc B litz (1993), and the values for the Galacto- 
centric distance (7.6±0.3 kp c) and Solar angular v elocit y 
(214±7km s"^) compiled by iKothes fc Doughertvl l|2007f ). 
It can be seen from the plot that the average G305 radial 
velocity is close to the tangential point. Once the uncertain- 
ties in the rotation curve are taken into account, we find a 
distance of 4.2±2.0 kpc. 

3.2.2 Spectrophotometric distances of the two clusters 

In calculating the spectrophotometric distances to each clus- 
ter, we begin by analysing those stars for which we are con- 
fident of the luminosity classes. As discussed in Sect. 12.2.31 
the star Dl-3 shows the clear spectral sig natures of a 08- 
B3 supergiant. From [Crowt her et al.l l|2006l ). such stars have 
absolute K-hand magnitudes of Mk = — 6.2±0.5, and so for 
this star we find a spectrophotometric distance to Danks 1 of 
3.48lQ'7jkpc. Applying the same analysis to the star D2-1, 
another 08-B3 supergiant, we find a distance to Danks 2 of 
2.6llQ gjkpc. Both these distances are within the errors of 
one another, as well as being consistent with the kinematic 
distance derived in the previous section. 

In Table [4] we list the spectrophotometric distances for 

^ |http : //www ■ ast ■ leeds ■ ac ■ uk/RMS | 



the remaining stars in each cluster which do not display 
evidence for high luminosity classes (i.e. no Br 7 emission, 
indicative of supergiants). Distances are calculated for each 
of the V, in and I luminosity classes, based on the absolute 
-R'-band ma gnitudes quoted i n the 'observational' tempera- 
ture scale of Martins fc Pled (|2006l ). By making the a priori 
assumption that the stars with similar spectral types have 
the same luminosity class we can assign luminosity classes 
that form a consistent picture of the distances to each clus- 
ter. 

In Danks 1, consistent distances are found for all the 
'non-supergiants' if we assign them to be class v objects, 
i.e. on or near the main sequence. A weighted mean of these 
distances gives 4.16 ± 0.6kpc, where the uncertainty is the 
r.m.s. deviation on the mean. Again this is consistent with 
that derived from the supergiant in this cluster and the kine- 
matic distance. 

In the case of Danks 2, the 06-8 stars all have Jf-band 
fluxes which are consistent with them being dwarfs, when 
taking into account both the kinematic distance and the 
spectrophotometric distance of D2-1. These stars have an av- 
erage spectrophotometric distance of 3.4±0.2kpc, consistent 
with G305's kinematic distance of 4.2±2.0kpc. The 08-9 
stars however appear to be too bright to be class v stars, and 
have probably evolved away from the main- sequence. There 
K-hand brightnesses are more typical of giants, though we 
consider the absolute brightnesses of class in stars to be too 
uncertain to determine a reliable spectrophotometric dis- 
tance. 

Since the distances to the two clusters are within the 
errors of one another, and the distance of Danks 2 is con- 
strained by only two measurements, for the rest of this paper 
we make the assumption that the two are at the same dis- 
tance. This is a reasonable assumption to make, since if the 
size of the G305 cloud along the line-of-sight is comparable 
to its angular size, the maximum difference in distances be- 
tween the two clusters is only ~30pc. Taking the weighted 
average of the spectrophotometric distances for the dwarfs 
only, as well as the kinematic distance, we find a distance to 
the clusters of 3.8±0.6kpc. 

3.3 Stellar populations and cluster ages 

The stellar populations of the two clusters - the presence 
of early-to-mid O dwarfs and supergiants - clearly indicate 
ages <6Myr. At these young ages it is not useful to simply fit 
model isochrones to the near-IR colour-magnitude diagrams, 
as there is a large degeneracy in age. Instead, we estimate 
ages for the two clusters by three methods: analysis of their 
stellar populations; examining the MS turn-offs; and also by 
studying the low-mass pre-MS population of each cluster. 

3.3.1 Danks 1 

In the previous section we argued that the 04-6 stars in 
Danks 1 were very likely on the main-sequence, whereas the 
stars which are slightly brighter (the 06-8If, WNLh and OBI 
stars) appear to be post-MS objects. From this information, 
we conclude that the MS turn-off corresponds to a spectral- 
type of 04-6. From Fig. O this implies a ZAMS mass of 
6O±2OM0 once the error in distance is taken into account 
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Table 4. Spectrophotometric distances of the stars in each cluster, assuming luminosity classe s V, III and I , deno ted as D^, Dm, 
respectively. Calculations use the absolute X-band mag nitudes of Martins & Plcz (2006) and Cr owther et alj l|2006l l. 



Star m222 Spec Type Dv/kpc Dm/kpc Di/kpc 



Danks 1 




















Dl-7 


9.511 


04-06 


3.21 




4.15 


4.59 


•H- 


5.20 


6.08 ^ 6.11 


Dl-8 


9.655 


04-06 


3.43 


■H- 


4.44 


4.91 




5.56 


6.50 •(-!> 6.53 


Dl-9 


10.212 


04-06 


4.43 


■H- 


5.73 


6.34 




7.19 


8.40 -H- 8.44 


Dl-10 


10.581 


06-08 


4.06 


•H- 


5.25 


6.58 




7.52 


9.96 ^^ 9.96 


Dl-11 


11.229 


09-B3 


1.19 


•H- 


4.78 








14.18 24.42 


Dl-12 


11.569 


09-B3 


1.39 


«• 


5.60 








16.59 ^ 28.56 


Danks 2 




















D2-4 


9.510 


08-09 


2.37 




2.70 


4.11 




4.39 


6.64 ^ 7.02 


D2-6 


9.606 


08-09 


2.47 


•H- 


2.83 


4.30 




4.58 


6.94 .H. 7.33 


D2-7 


9.644 


08-09 


2.52 


•H- 


2.88 


4.37 




4.67 


7.06 ^ 7.46 


D2-5 


9.583 


06-08 


2.80 


■(->■ 


3.62 


4.53 




5.18 


6.86 6.87 


D2-8 


9.805 


06-08 


3.10 


■H- 


4.01 


5.02 




5.74 


7.60 ^ 7.60 



ijMevnet fc Maedeijl2000l : lMartins et al.ll2005ll!|. The MS life- 
time for such stars is around 3-4Myr ijMevnet fc Maederl 
I2OOOI ). which is therefore an upper limit to the age of 
Danks 1. 

In the CMD of Fig. \E\ a departure from the ZAMS is 
seen at magnitudes fainter than F222M~15. This can be 
interpreted as the point at which lower mass stars in the 
cluster are beginning to arrive on the MS. From the implied 
mass at which stars are on the MS, we can get an indepen- 
dent measure of the cluster's age. The brightness at which 
the stars are seen to join the MS corresponds to a mass of 
4±l Mp). From comparison with models of p re-MS evolution 
(e.g. lPalla fc Stahleij[T999l : ISiess et al.ll20od ). we find the age 
of the low-mass stellar population of Danks 1 to be l-2Myr. 

The most massive stars present in the cluster can also 
be used as an age discriminant, since the lifetime of a star is 
closely related to its initial mass. The most massive stars in 
Danks 1 are likely to be the three WNLh stars. Such stars are 
also present in the Arches cluster, where their initial masses 
have been found to be as high as I2OM0 (JVlartins et al. 
|2008| ). Following Martins et al., we use their mean A"-band 
bolometric correction for WNLh stars, BCk = —4.21 ±0.26, 
as well as the average cluster extinction fSect. [371]) and dis- 
tance modulus to G305 (Sect. 13. 2p . Using these numbers, 
we find that the luminosities of Dl-1, Dl-2 and Dl-5 are 
log(I//L0) = 6.5, 6.2, and 5.9 respectively (all ±0.2dex). 
Clearly these three stars are all intrinsically luminous, and 
therefore very massive. While accurate spectrophotometric 
mass determinations await quantitative modelling of their 
spectra, we can tentatively say here that the masses of Dl-1 
and Dl-2 are likely to be in excess of 90Mb (com paring with 
Geneva rotating models. [Mevnet fc Maedeill200Cll ). This then 
places an upper limit to the age of the cluster of 3Myr. 

To summarize, an age of for Danks 1 is con- 

sistent with all pieces of evidence from both the high mass 
stars and the low mass pre-MS stars. While it has been found 



' We determine a scale of spectral-type versus mass using the 
mass-temp erature scale defined by the stellar structure models of 
llMevnet fc Ma cdcr 2000), and the spectral- type - temperature 
scale of lMarUns et al.l | |2005|) . 



that these two age indicators can give contradictory results 
(e.g. Westerlund l,,Brandner et al. 2008;) . and that pre-MS 
isochrones may give ages which a re systematic ally younger 
than those indicated by MS stars l|Navloill2009l 'l. we find no 
evidence for such a discrepancy here. 

Finally, we mention a star in Danks 1 that does not seem 
to fit with our derived age for the cluster. The star Dl-3 has 
the appearance of a 'normal' blue supergiant, with narrow 
absorption lines, and the absence of significant line emission 
implies a relatively weak wind. It is therefore natural to 
conclude that this star is one of moderate initial mass, say 
~2O-4OM0, in an advanced evolutionary state. This would 
imply an age of 4-lOMyr, which is clearly at odds with the 
other evidence from the cluster. This suggests that Dl-2 was 
not born with the rest of the stars in Danks 1, and is instead 
part of a population of older stars which includes the other 
evolved massive stars seen in the field of G305 (see Table [3]). 
It is possible that these stars formed along with the older 
Danks 2 (see next section) and were dynamically ejected. 
For a projected cluster separation of 40pc, if the star was 
ejected from Danks 2 ~3Myr ago then this implies a runaway 
velocity of ~15km s~^, which is certainly not unreasonable. 

3.3.2 Danks 2 

Using the same diagnostics of the cluster age, the evidence 
suggests that Danks 2 is somewhat older than Danks 1. The 
earliest spectral type MS stars are 06-8, consistent with 
masses of 3O-4OM0 and hence an upper limit to the age of 
4-6Myr. The point at which the pre-MS stars join the MS 
is at a lower mass than in Danks 1 - ~2M0, indicating an 
age of 3Myr with upper and lower limits of 2-lOMyr due to 
the uncertainty in distance. 

The most luminous stars in Danks 2 can again be used 
as an age indicator, though the results are less conclusive 
than in Danks 1. In particular, carbon-rich WR stars should 
be present in clusters with ages between 3-6Myr, according 
to rotating Geneva models (Mcynct & Maedcr 2000). The 
yellow star, D2-2, lies to the left of the ZAMS track in the 
colour- magnitude diagram of Fig. [51 suggestive that the star 
does not belong to the cluster. 

Our inability to determine a precise spectral type for 
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the OB supergiant D2-1, the brightest star in Danks 2 in 
the near-IR, means that it places only weak constraints 
on the cluster age. The if -band bolometric correction of 
a 08- B3 supergiant is BC k = — 3.93 —1.55, accord- 
ing to iMartins fc Pie3 (|2006l ) and ICrowther et~aLl (|2006l ). 
which implies a bolometric luminosity of log(LD2-i/io) ~ 

5.3 — >■ 6.3 for the coolest and hottest temperatures respec- 
tively. A luminosity for D2-1 closer to the lower limit seems 
more likely, since more luminous stars do not tend to have 
the spectroscopic appearance of 'normal' blue supergiants, 
and instead have strong emission lines owing to their dense 
winds. 

From the above evidence, we suggest that the age of 
Danks 2 is 3t?Myr. The error is dominated by the uncer- 
tainty on the cluster's heliocentric distance. However, we 
can be certain that the two clusters are at the same dis- 
tance, within a few percent. Therefore, it seems very likely 
that Danks 2 is the older of the two clusters, being created 
1.5lJ;5Myr before Danks 1. 

3.4 Initial mass functions 

In order to construct initial mass functions (IMFs) for the 
two clusters, we must first create luminosity functions (LPs) 
which are decontaminated of foreground stars. We experi- 
mented with two methods of doing this. In the first method, 
we contructed LPs for both the cluster and control fields, 
then subtracted one from the other once completeness ef- 
fects had been taken into account. Secondly, we used the 
P222M magnitudes of the stars in the field-corrected colour- 
magnitude diagrams of Pigs. [5]and[6l The completeness cor- 
rections applied to these LPs were those of the cluster fields, 
as these completeness limits are typically higher than those 
of the control fields. In any event, the limit down to which we 
chose to measure the IMPs is well above the point at which 
the photometry becomes incomplete (see below). We also 
experimented with de-reddening the stars in the decontam- 
inated LP back onto the ZAMS track, under the assump- 
tion that stars lying slightly to the right of the track are 
cluster stars with extra local extinction. Any stars which 
had (P160W-P222M) excesses greater than O.Bmags were 
discarded, as these most likely belong to a separate back- 
ground population, though we note that without spectra we 
cannot completely rule out that these objects are local to 
the clusters. 

To convert the LPs to IMPs, we use the evolutionary 
models of the Geneva group. The lates t versions of these 
take into account stellar rotation (e.g. iMevnet fc Maedeil 
l200dl ). but computations do not exist for stars with masses 
below 9M0, since rotation does not have a large impact on 
intermediate to low mass stars. Por this reason we splined 
together the rotating models with the older non-rotating 
models of dSchaller et al.iil992 ) at 9M0. As a check on the 
robustness of our results, we also experimented with the 
older non-rotating models with higher/lower metallicities, 
and varying mass-loss rates. Instrumental magnitudes and 
colours were determined the same way as in Sect. 13.31 

In Pig. [8] we show the IMPs for the two clusters, de- 
rived using the colour-magnitude corrected LPs, and using 
Solar-metallicity isochrones including rotation at the ages 
indicated in each panel. In each case, we measure the slope 
of the IMP only in the mass ranges where we are sure that 



stars are on the main-sequence. The low mass stars which 
have not yet reached the MS will have colours and magni- 
tudes which are systematically different from those predicted 
by main-sequence evolution codes; while the post-MS be- 
haviour of high-mass stars is extremely uncertain. The pre- 
MS mass limit is measured from the kink in the CMD seen in 
both cluster sequences at mp222 ~ 15. The post-MS mass is 
determined from the point at which the stellar spectropho- 
tometric distances are no longer consistent with luminosity 
class V stars. In practise, the upper limit for the range in 
which we measure the IMP (dotted blue lines in Pig. [8]) is 
set by the point at which the number of stars per bin drops 
below 2, where the uncertainties become non-poissonian. 

In each cluster we see slopes that are consistent with 
the Salpeter value of F = —1.35. Danks 1 has a slope 
F = -1.40±0.17, while Danks 2 has a slope F = -1.38±0.16. 
These measurements are robust to the type of evolutionary 
model used in the isochrones; the effect of varying mass- 
loss rate, stellar rotation and metallicity is small compared 
to the uncertainty. Varying the cluster age within the mea- 
sured uncertainties has little impact, since we measure the 
slope of the IMP only for objects on the MS where stellar 
properties change very little over ~10''yrs. Finally, varying 
the cluster distances between the upper and lower limits 
(3.2-4.4kpc) produces changes in the slope which are small 
compared to the uncertainties. We can therefore say that we 
find no significant evidence for variations in IMP between 
the two clusters in G305, and no significant deviation from 
the Salpeter slope. 

If, as we have argued in Sect. 13.31 Danks 1 is in a pre- 
supernova (SN) state, the IMF of this cluster may be used 
to investigate the upper end of the IMF. Observations of the 
Arches cluster have indi cated that t here may be an upper 
mass limit of ~150Mo (jPigeil 120051 ). though observations 
of R136 in the Large Magellanic Cloud su ggest that stars 
may form with masses greater than this (jCrowther et al.l 
I2OI1J ). Quantitative spectral modelling of the brightest stars 
in Danks 1 would yield accurate bolometric luminosities, 
and, by comparing to evolutionary models, initial masses. 
It would then be possible to discuss these objects in the 
context of the high-mass end of the IMP. 

3.5 Cluster masses 

In order to determine the initial masses of the two clusters, 
we simply fit the observed IMPs measured in the previous 
section with a functional form. Since all functional forms 
of the IMP have a slope which is approximately Salpeter- 
like in the mass-range appropriate to our measurements, we 
simply scale the functional IMP to fit our observations and 

then integrate over all stellar masses. 

Fitting a Kroupa IMP to our observations (jKroupal 

we find that the initial masses of Danks 1 and Danks 2 
are 8OOO±15OOM0 and 3OOO±8OOM0 respectively. The un- 
certainties are due to a combination of errors carried forward 
from cluster age and distance, as well as a stochastic error 
due to low number statistics at the high-mass end. 

We note that the Kroupa IMP puts a large amount of 
mass into stars with sub-Solar masses, going all the way 
down to O.OIM0. If we were to use, for example, a Chabrier 
IMF (|Chabrierll200l ). we would find that the inferred cluster 
masses would be approximately a factor of two lower. 
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Figure 8. Initial Mass Functions (IMFs) for the two clusters: Danks 1 {left) and Danks 2 {right). The solid and dotted histograms 
show the data before and after correction for completeness. The green dashed line illustrates the Salpeter slope. The slope of each mass 
function, shown as the red solid line, is computed between the MS turn-on and turn-off marked by the blue dotted lines. 



Table 5. Summary of the clusters' physical properties 



Cluster 


Age (Myr) 


Minit (M0) 


R0.5 (pc) 


log(p/M0pc-3) 


Danks 1 
Danks 2 


^■°-0.5 

OQ + 3.0 


8000±1500 
3000±800 


0.17±0.05 
0.36±0.09 


5 .5+0.5 
4 2+0.5 
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Figure 9. Radial surface brightness profiles for the two clusters. 



Fig. [9l The error bars on each point represent the effect 
of varying the position defined as the cluster centre. The 
half-hght radii of the clusters are given in Table [3 where 
the errors on these values take into account both the error 
on i?o.5 and the uncertainty in distance. Also summarized 
in Table [5] are the cluster densities, ages and intergrated 
masses. 



4 DISCUSSION 

Our photometric and spectroscopic analysis of Danks 1 & 2 
reveal that they are both rather massive, young and com- 
pact, with a total mass of > lO^M©, and both forming over a 
short interval of ~2 Myr. In the following section we discuss 
how these clusters - and by extention the G305 complex - 
compare to other Galactic and extra-galactic star-forming 
regions. 



3.6 Cluster sizes and densities 

We determined the half-light radius of each cluster R0.5 by 
measuring the cumulative surface brightness profiles (SBPs) 
as a function of distance from the cluster centre. Since our 
NICMOS observations have a field-of-view too small to de- 
termine at which radius the SBPs fall to the ambient level, 
we used 2MASS images to make these measurements. We de- 
fined the cluster centre as the centre of the NICMOS fields- 
of-view, but since the cluster centre is not easily defined we 
repeated the measurements with the centre offset in RA and 
DEC by ±3". 

The cumulative SBPs of the two clusters are shown in 



4.1 Cluster comparison 

With a mass of ~8000 Mq, an age of ~1.5 Myr and a 
stellar population containing both WNLh and early-mid 
0-type supergiants, Danks 1 closely resembles a number 
of other massive young clusters which are also associated 
with Giant Hll regions. These in clude Trump ler 14 (lo- 
cated withi n the Carina nebula; ISmith et al.l 2006) and 
NGC 3 603 jHaravama et al]|2008h . as weU as the Arches 
cluster (|Figer et al.ll2002l l. which, given its location within 
the central 50 pc of the Galaxy, cannot be uniquely associ- 
ated with a natal birthcloud. 
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By comparison, Danks 2 is both older, less massive and 
less dense than the above clusters. Its mass and stellar pop- 
ulation (and hence age) are reminisce nt of NGC 6231 at the 
centre of Sco QBl (fSung et al.|[l993 ). or CI 1806-20 (|Figeil 
I2OO5I : iBibbv et al.ll2008l '). though the latter cluster does not 
have any obvious associated GMC. 

The difference in densities between Danks 1 and 
Danks 2 may be related to the dynamical evolution of young 
massive star clusters. The Central Cluster aside, Danks 1, 
the Arches, and NGC 3603 all have ages which are thought 
to be below 3Myr, and so all may be yet to experience a SN 
of a cluster member. The removal of mass from the cluster 
by SN e may leave the cluster su per- virial, causing it to ex- 
pand (|Hillslil980l : lGoodwin fc Bastian.2006 '). If Danks 2 has 
already experienced SNe, this may explain the difference in 
central density between the two clusters. 



4.2 Cluster complex comparison 

If Danks 1 closely resembles both Trumpler 14 and 
NGC3603, one might also ask whether the properties of the 
associated star forming regions are also comparable? At this 
point it is instructive to extend this comparison to 30 Do- 
radus and its central cluster R136 in the Large Magellanic 
Cloud, the most luminous Giant H 11 region in the Local 
Group. 



4.2.1 Multiple stellar populations 

The evidence presented here suggests that star formation 
activity within the G305 commenced within the last 6 Myr. 
Following identica l arguments to tho se advanced for the 
W51 complex by IClark et all |2009a'), the apparent lack 
of luminous Red Supergiants (RSGs) within the complex - 
which would clearly be visible at IR wavelengths - is consis- 
tent with the picture that star formation was absent beyond 
~6 Myr ago. 

At least two further populations of massive (post- 
)MS stars are present within G305, with significant evi- 
dence for mo r e rece nt and ongoing massive star formation. 
iLeistra et al] (|2005l ) find a mid-O supergiant and two late- 
O/early-B main sequence objects located within the cluster 
they designate as G305-f00.2. Found within a small bubble 
on the periphery of the complex, Hindson et al. (in prep.) 
suggest an age of <1 Myr for the compact Hll regions asso- 
ciated with it. 

Finally, a substantial diffuse population of ma ssive 
stars appears to be p resent, with IShara et al.l (|2009l ') and 
iMauerhan et al.l (|2011h discovering additional WRs within 
the confines of the complex in addition to the members 
of Danks 1 & 2, WR 48a and the candidate dusty WCL 
MSX305. 4013+00.0170 (Table [3]). Indeed this WR popu- 
lation is significantly larger than that found within both 
Danks 1 & 2 combined, while the presence of four WC stars 
implies a minimum age of at least 3 Myr for this diffuse 
stellar component. 

The NGC3603 and 30 Dor and Carina star forming re- 
gions also all appear to host multiple stellar populations. 
Dominated by the young massive cluster NGC 3603 (1- 
2 Myr), the presence of the blue super giants Sher 23 fc 2 5 
reveals an older (~4 Myr) population i Melena et al.ll2008f ). 



(IWalborn fc Blades|[l997l ) likewise report a number of differ- 
ent stellar groups within 30 Dor; (i) the central cluster R136 
(2-3 Myr), (ii) an older (~4-6 Myr) population distributed 
across the region, (iii) Hodge 301; a ~10 Myr cluster to 
the NW of R136 (iv) the R143 association in the south-east 
(~6-7 Myr) and (v) a very young (~1 Myr) population on 
the pe riphery of the com plex, which we return to below. 
Lastly, 'S mith et all (|2006l ') summarise the properties of Ca- 
rina, again reporting a range of ages for the clusters located 
within it (Trl4 at ~1-1.5 Myr, Trl5 at 6±3 My r, Trl6 at 2- 
3 Myr, BolO at ~7 Myr and Boll at <3 MvrV lSmith et all 
(2006) also comments on the fact that in contrast to the 
previous examples, these populations are rather widely dis- 
tributed over the complex. While the richest clusters (Trl4 
& 16) are centrally located the remainder are observed at 
significant distances (>10 pc)in the outer regions of the com- 
plex. An age spread amongst the st ellar populations w ithin 
the W51 complex is also apparent (|Clark et al.ll2009al ). al- 
though in this case there is no evidence for a compact central 
cluster, with star formation apparently distributed through- 
out the host GMC. 



4.2.2 Ongoing star formation in cluster complexes 

As reported by IClark fc Ported l|2004l ). IClark et all (|201lD 
and iHindson et all (|2oioi r there is compelling evidence for 
a further generation of massive stars forming with the G305 
complex. Methanol masers and ultracompact (UC-)IIn re- 
gions, both unambiguous indicators of ongoing massive star 
formation, are present within G305, being predominantly lo- 
cated on the periphery of the bubble. Integrated IR- and 
radio-fluxes likewise argue for the presence of a signifi- 
cant population of embedded massive YSOs, while sub-mm 
observations reveal the presence of a substantial reservoir 
(~6xl0^ Mq) of cold molecular gas available to fuel future 
activity. 

The star forming complex assoc iated with NGC 3603 
also s hares these properties, with iNiirnberger fc Stankd 
^200^ describing the presence of methanol massers, dusty 
embedded sources and molecular cores within the confines 
of the remnant natal GMC. Both IR- and radio-fluxes are 
likewise com parable to those of G305 to within a factor of 
a '^few (e.g. ICrowther fc Contill2003l ). A comparable mor- 
phology is also observed for 30 Dor, with active star for- 
mation loc ated on the periphery o f the cavity surround - 
ing R136 (|Walborn fc Blades! 1 19971 : IWalborn et all I2OO2I '). 
Finally, while such spatial segregation between pre- and 
(post-)Main Sequence stars is less apparent with in Carina, 
the presence of ucHll regions (| Brooks et al.|[200ll ) points to 
onoing massive stars formation, pa rticularly in the 'South - 
ern Pillars' region of the complex (|Rathborne et al] 120041 '). 
Moreover, both integrated IR- fc radio-fluxes are directly 
comparable to those of G305, as is t he mass of cold molec- 
ular material (|Smith fc Brooks! [2OO7I ). 

In all four cases it has been suggested that the com- 
plex morphologies reflect the propogation of triggered star 
formation through the molecular cloud. This impression is 
particularly strong for G305, where the active star form- 
ing regions appear restricted to the periphery of the cavity, 
and indeed with a 'third generation' of masers and ucHll 
regions associat ed with the bubble surrounding the clus - 
ter G305-I-00.2 l|Clark fc Portei!l2004l : iHindson et al.ll2O10! '). 
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While such morphologies are necessary, they are not con- 
sidered sufficient evidence for triggered star formation since 
it might be supposed that the radiation- or wind blown- 
bubble is simply uncovering exisiting activity. In this re- 
gard our unbiased multiw avelength datasets , encompassing 
the whole G305 complex (|Clark et al.ll201ll ) will be partic- 
ularly valuable in distinguishing between these possibilities 
via the identification of (massive) YSOs. If star formation 
was sequential, one would expect to find YSOs interior to the 
cavity but none outside, whereas they should b e distributed 
throughout the GMC if this was not the case fcf. lSmith et alj 
Indeed, our dataset will allow us to constrain not only 
the propagation of star formation through the GMC, but 
also the relative stellar yields of each successiv e generation, 
a par ticularly interesting prospect given that ISmith et alJ 
l|2010 l) suggest that this decreases with time (and successive 
generations) in Carina. 

Despite having a comparable cloud mass (Parsons et 
al., accepted with minor revision) and massive star for- 
mation having been underway for several million years 
(cf. radio- and I R-fiuxes similar to the above regions; 
IConti fc Crowtheil i2004'). it is interesting that the W51 com- 
plex appears not to host such triggered star formation. In- 
stead, multiple, causally unconnected, star forming regions 
are found throughout the cloud, albeit possibly synchro- 
nised by the presence of an external ag ent such as the pas- 
sage of a Galactic spiral density wave (| Clark et al. I l2009al : 
iKumar et al.|[20o3 ). Nevertheless, the physical endpoint in 
terms of the production of multiple young massive clusters 
with an age range of < 10^ years appears remarkably similar 
to the above complexes. 



4.3 A template for star formation 

Irrespective of the mechanism of star formation (multiseeded 
versus triggered) , these complexes present a template of how 
a IO^-IO^Mq GMC is converted into stars. The process is 
distributed in both time and space, yielding a number of 
rather massive (1O"^-1O*M0) clusters with a significant age 
spread over a region spanning an angular diameter of ~30pc. 
Of immediate interest is the fact that the properties of both 
NGC 3603 and Danks 1, in terms of integrated mass, den- 
sity, age and stellar content is directly comparable to that 
of the Arches. Clearly the production of such clusters is not 
dependent on the extreme conditions of the Galactic centre, 
although with 3 clusters of similar masses (a few xlO^M©; 
the Arches, Quintuplet and Galactic Centre) the star for- 
mation rate here appears to have been larger over the past 
few Myr than in the NGC 3603 and G305 complexes. 

Globally, these giant star-forming regions are also of 
interest due to the fact that they mirror the structure of cx- 
tragalactic star cluster complexes, albeit with cluster masses 
several orders of magnitude smaller. Indeed the integrated 
spectral energy distribution of Carina is rem arkably sim- 
ilar to those of U l tralum inous IR Galaxies (|Smith l2007l : 
ISanders fc Mirabej [199^ ') . suggesting that a similar mode 
of star formation is also present in such galaxies. Because 
of this, determining the energy budget of Galactic examples 
such as Carina and G305 is important since it will permit 
a calibration of the star formation rate from the emergent 
fiux. 

While an accurate determination of the energy budget 



of G305 requires both accurate stellar parameters and the 
mid-far IR fiux (Clark et al. in prep.), we simply note here 
that the current feedback from Danks 1 fc 2 is likely to 
be dominated by the 3 WNLh stars. Utilising the Lyman 
flux estimates for similar WNLh stars in th e Arches clus- 
ter sc aled to the luminosities of these objects (|Martins et al.l 
|200j) , the amount of ionizing radiation emitted by the three 
WNLh stars in Danks 1 alone appears comparable to the to- 
tal ionizing fiux within the entire G305 region estimated by 
I Clark fc Porteil (|2004l '). This suggests that, as with Carina, 
G305 suffers significant photon leakage. Hence, estimating 
the stellar contents and star formation rates of such regions 
using either IR or radio fiuxes alone may result in signifi- 
cant underestimates for all but the very youngest complexes, 
where feedback has yet to uncover embedded clusters. 



4.4 Future evolution 

Finally we turn to the long term evolution of the G305 com- 
plex. As stellar evolution drives mass loss, the long term sur- 
vivability of the clus ters (Danks 1 fc 2 and G3 05-I-00.2) is 
somewhat uncertain ijCoodwin fc Bastianll200^ 'l. If they do 
disperse rapidly, G305 will increasingly resemble a classical 
OB association, noting that the stars within such aggregates 
appear to have formed over a comparable period (e.g. Cyg 
OB2 fc Sco OBI; Clark et al. submitted: iNegueruela et aU 
|2008| ). However, at present Danks 1 fc 2 appear to be tightly 
bound, with ratios of ages to crossing times of order ~10, 
suggesting tha t they will survive a s clusters for several mil- 
lion years (see lGieles fc Portegies Zwart,2011, V In this case, 
an obvious point of comparison is the Perseus complex, 
which contains two distinct star clusters, h and x Persei, 
surrounded by a halo of stars of simil ar age. The comple x 
has been recently studied in detail bv ICurrie et all (|2010f ). 
finding remarkably similar ages and masses for the two clus- 
ters, 3700-4200MGlf| and 14±lMyr, and an age for the halo 
population which was only marginally younger, 13.5±lMyr. 

A quantitative comparison with our results presented 
here show that this is somewhat different to what we see 
in G305. The age spread we find for Danks 1 and Danks 2 
is small, 1-2 Myr, and will become more difficult to detect 
as the clusters approach the same age as h and x Persei. 
However, we would still expect to be able to detect a sig- 
nificant age difference between the central clusters and the 
surrounding halo of stars. It is clear that the remains of the 
G305 cloud is still forming stars today, implying an age dif- 
ference of up to 5 Myr. Therefore we may still expect to 
detect a difference in the ages of the clusters and the sur- 
rounding halo when the clusters are well beyond 20 Myr 
old. 

An alternative comparison might be with the associ- 
ation of RSG dominated clusters at the base of the Scu- 
tum Crux arm, which appear to have f ormed over an ex- 
tended period betw een ~ 12-20 Myr ago (jPavies et al. I I2OO8I : 
IClark et al.ll2009bl ). In particular, RSGC3 would appear to 
form a particularly good comparator, being surrounded by 
at least one further lower mass cluster and a number of 
isolated RSGs of similar ages and masses within a ~30 pc 



* These authors assumed a Miller-Scalo IMF to determine the 
total masses of the populations in the Perseus complex. 
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radius (|Clark et al.ll2009bl : iNegueruela et al.ll2011bl ). While 
currently separate, with a (projected) separation of only 
~4pc a combination of cluster expansion driven by stel- 
lar mass loss and dynamical interaction could lead to the 
effective merger of Danks 1 & 2, further emphasising the 
similarity between the 2 regions. 

In the case of Danks 1& 2 if this were to occur before the 
RSG dominated phase (>10 Myr) this would lead to a mix- 
ture of stellar spectral types not predicted by current evolu- 
tionary theory. A similar explanation has been proposed to 
explain the stellar population of the young (~7Myr) massive 
cluster NGC 1569-A.The integrated spectrum displays fea- 
tures attributed both to Red Supergiants (RSGs) a nd WRs 
iGonzalez Delgado et all Il997l : iHunter et al] I2000I '). which 
stellar population models struggle to produce simultane- 
ously for a single instantaneous starburst. Analysis of the 
point-spread function and velocity dispersion has also indi- 
cated two d istinct components to NGC 1569-A, now labelled 
Al and A2 (Ide Marchi et al.|[l997l : [Gilbert fc~Graham. 20021 
Further. iMaoz et al.l ( 2001 ) found that the WR spectra sig- 
natures were confined to A2. The current understanding of 
NGC 1569-Al and NGC 1569-A2 is that they have compa- 
rable masses of ~ 6 x 10^ M©, but ages of ;^5Myr and ~8Myr 
respectively; a comparable age spread to that found between 
Danks 1 and Danks 2, though the NGC 1569-A clusters are 
slightly older and orders of magnitude greater in mass. 



5 SUMMARY AND FUTURE WORK 

We have provided the first comprehensive study of the two 
central clusters of the G305 star-forming complex, Danks 1 
and Danks 2. We determine a distance to the clusters and 
the host molecular cloud of 3.8±0.6kpc. This takes into ac- 
count both the spectrophotometric and kinematic measur- 
ments, which are in agreement with each other. The total 
stellar masses of Danks 1 and Danks 2 are SOOOilSOOM© 
and 3OOO±8OOM0 respectively. Analysis of their initial 
mass functions (IMFs) shows that both are consistent with 
Salpeter. A further population of apparently isolated mas- 
sive stars, predominantly comprising WRs, is also present 
within the confines of the complex. 

From analysis of the stellar content and pre-main- 
sequence stars of the two clusters, we have been able to 
piece together the star-forming history of the G305 com- 
plex. Danks 2 is the oldest object in the complex, with an 
age of 3l"jMyr. The younger age of Danks 1, l.SlJ ^Myr, is 
consistent with its formation being triggered by the feedback 
of Danks 2. It is likely that the combined winds from the two 
clusters are responsible for the evacuation of the complex's 
central cavity, and perhaps are responsible for the consid- 
erable amount of star formation that has occurred in the 
last ~0.5Myr around the periphery of the G305 cloud. The 
origin of the diffuse population of evolved massive stars is 
less clear, though the presence of both WN and WC stars is 
consistent with this population forming at a similar time to 
the two central clusters. 

The masses and stellar content of the two clusters, com- 
bined with the significant numbers of Wolf-Rayet stars in 
close vicinity to the clusters, and the numerous massive 
protostars in the surrounding molecular cloud, make G305 
one of the most bountiful regions for massive star forma- 



tion known in the Galaxy, comparable to the Carina nebula. 
As such, the region represents the perfect template for cali- 
brating unresolved observations of similar regions at extra- 
galactic distances. 

In terms of work for the future, the high-quality spec- 
tra presented here will be modelled with stellar atmosphere 
codes. This will allow a qualitative estimate of the amount of 
feedback from the clusters, both in terms of the ionizing ra- 
diation and the mechanical energy from stellar winds. Also, 
it will provide an estimate of the bolometric luminosities for 
the most massive stars in the cluster, and the universality 
of the IMF up to the highest stellar masses. 



ACKNOWLEDGMENTS 

We thank Fabrice Martins for supplying us with the spectra 
of stars in the Arches cluster. We also thank the anonymous 
referee whose comments and suggestions helped improve the 
paper. BD is supported by a fellowship from the Royal As- 
tronomical Society. This work is in part based on observa- 
tions made with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. These 
observations are associated with program #11545. Support 
for program #11545 was provided by NASA through a grant 
from the Space Telescope Science Institute, which is op- 
erated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. This 
work is in part based on observations collected at the Euro- 
pean Organisation for Astronomical Research in the South- 
ern Hemisphere, Chile, under programme number 077. C- 
0207(B). Financial support from the Spanish Ministerio de 
Ciencia e Innovacion under projects AYA2008-06166-C03-02 
and AYA2010-21697-C05-01 is acknowledged. 



REFERENCES 

Alves, J. & Homeier, N. 2003, ApJ, 589, L45 

Bastian, N., Emsellem, E., Kissler-Patig, M., & Maraston, 

C. 2006, A&A, 445, 471 
Bastian, N., Gieles, M., Efremov, Y. N., & Lamers, 

H. J. G. L. M. 2005, A&A, 443, 79 
Bibby, J. L., Crowther, P. A., Furness, J. P., & Clark, J. S. 

2008, MNRAS, 386, L23 
Bohannan, B. & Crowther, P. A. 1999, ApJ, 511, 374 
Bonanos, A. Z., Stanek, K. Z., Udalski, A., Wyrzykowski, 

L., Zebruh, K., Kubiak, M., Szymariski, M. K., Szewczyk, 

O., Pietrzyhski, G., & Soszyhski, I. 2004, ApJ, 611, L33 
Brand, J. & Blitz, L. 1993, A&A, 275, 67 
Brandner, W., Clark, J. S., Stolte, A., Waters, R., Negueru- 

ela, I., & Goodwin, S. P. 2008, A&A, 478, 137 
Brooks, K. J., Storey, J. W. V., & Whiteoak, J. B. 2001, 

MNRAS, 327, 46 
Chabrier, G. 2003, PASP, 115, 763 

Clark, J. S., Davies, B., Najarro, F., MacKenty, J., 
Crowther, P. A., Messineo, M., & Thompson, M. A. 2009a, 
A&A, 504, 429 

Clark, J. S., Davies, B., & Thompson, M. A. 2011, Bulletin 
de la Societe Royale des Sciences de Liege, 80, 223 



The young star- clusters Banks 1 and Banks 2 17 



Clark, J. S., Negucrucla, I., Davics, B., Larionov, V. M., 
Ritchie, B. W., Figer, D. F., Messineo, M., Crowther, 
P. A., & Arkharov, A. A. 2009b, A&A, 498, 109 
Clark, J. S. & Porter, J. M. 2004, A&A, 427, 839 
Conti, P. S. & Crowther, P. A. 2004, MNRAS, 355, 899 
Crowther, P. A. & Conti, P. S. 2003, MNRAS, 343, 143 
Crowther, P. A., Lennon, D. J., & Walborn, N. R. 2006, 

A&A, 446, 279 
Crowther, P. A., Schnurr, O., Hirschi, R., Yusof, N., Parker, 
R. J., Goodwin, S. P., & Kassim, H. A. 2010, MNRAS, 
408, 731 

Currie, T., Hernandez, J., Irwin, J., Kenyon, S. J., Tokarz, 

S., Balog, Z., Bragg, A., Berlind, P., & Calkins, M. 2010, 

ApJS, 186, 191 
Davies, B., Figcr, D. F., Kudritzki, R.-P., MacKenty, J., 

Najarro, F., & Herrero, A. 2007, ApJ, 671, 781 
Davies, B., Figer, D. F., Law, C. J., Kudritzki, R.-P., Nar 

jarro, F., Herrero, A., & MacKenty, J. W. 2008, ApJ, 676, 

1016 

Davies, B., Lumsdcn, S. L., Hoarc, M. G., Oudmaijer, 
R. D., & de Wit, W. 2010, MNRAS, 402, 1504 

de Marchi, G., Clampin, M., Greggio, L., Leitherer, C, 
Nota, A., & Tosi, M. 1997, ApJ, 479, L27+ 

Diolaiti, E., Bendinelli, O., Bonaccini, D., Close, L., Currie, 

D. , & Parmeggiani, G. 2000, A&AS, 147, 335 
Figer, D. F. 2005, Nature, 434, 192 

Figer, D. F., MacKenty, J. W., Robberto, M., Smith, K., 
Najarro, F., Kudritzki, R. P., & Herrero, A. 2006, ApJ, 
643, 1166 

Figer, D. F., McLean, L S., & Najarro, F. 1997, ApJ, 486, 
420 

Figcr, D. F., Najarro, F., Gilmore, D., Morris, M., Kim, 
S. S., Scrabyn, E., McLean, I. S., Gilbert, A. M., Graham, 
J. R., Larkin, J. E., Levenson, N. A., & Teplitz, H. L 2002, 
ApJ, 581, 258 

Gicles, M. & Portcgics Zwart, S. F. 2011, MNRAS, 410, L6 
Gilbert, A. M. & Graham, J. R. 2002, in lAU Symposium, 
Vol. 207, Extragalactic Star Clusters, ed. D. P. Geisler, 

E. K. Grebel, & D. Minniti, 471-+ 

Gonzalez Delgado, R. M., Leitherer, C, Heckman, T., & 

Cervino, M. 1997, ApJ, 483, 705 
Goodwin, S. P. & Bastian, N. 2006, MNRAS, 373, 752 
Hanson, M. M., Conti, P. S., & Rieke, M. J. 1996, ApJS, 

107, 281 

Hanson, M. M., Kudritzki, R.-P., Kenworthy, M. A., Puis, 
J., & Tokunaga, A. T. 2005, ApJS, 161, 154 

Harayama, Y., Eisenhauer, F., & Martins, F. 2008, ApJ, 
675, 1319 

Hills, J. G. 1980, ApJ, 235, 986 

Hiiidsoii, L., Thompson, M. A., Urquhart, J. S., Clark, 
J. S., & Davies, B. 2010, MNRAS, 408, 1438 

Hoare, M. G., Lumsden, S. L., Oudmaijer, R. D., Urquhart, 
J. S., Busfield, A. L., Sheret, T. L., Clarke, A. J., Moore, 
T. J. T., AUsopp, J., Burton, M. G., PurccU, C. R., Jiang, 
Z., & Wang, M. 2005, in lAU Symposium, Vol. 227, Mas- 
sive Star Birth: A Crossroads of Astrophysics, cd. R. Ce- 
saroni, M. FeUi, E. Churchwell, & M. Walmsley, 370-375 

Hunter, D. A., O'Connell, R. W., Gallagher, J. S., & 
Smcckcr-Hane, T. A. 2000, AJ, 120, 2383 

Johnson, H. L. 1966, ARA&A, 4, 193 

Koornneef, J. 1983, A&A, 128, 84 

Kothes, R. & Dougherty, S. M. 2007, ArXiv e-prints, 



0704.3073 
Kroupa, P. 2001, MNRAS, 322, 231 

Kumar, M. S. N., Kamath, U. S., & Davis, C. J. 2004, 

MNRAS, 353, 1025 
Larsen, S. S. 2004, A&A, 416, 537 

Leistra, A., Cotera, A. S., Liebert, J., & Burton, M. 2005, 
AJ, 130, 1719 

Maoz, D., Ho, L. C, & Sternberg, A. 2001, ApJ, 554, L139 
Martins, F., Hillier, D. J., Paumard, T., Eisenhauer, F., 

Ott, T., & Gcnzcl, R. 2008, A&A, 478, 219 
Martins, F. & Plez, B. 2006, A&A, 457, 637 
Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 

1049 

Mauerhan, J. C, van Dyk, S. D., & Morris, P. W. 2009, 

PASP, 121, 591 
Mauerhan, J. C, Van Dyk, S. D., & Morris, P. W. 2011, 

ArXiv e-prints 

Melena, N. W., Massey, P., Morrell, N. L, & Zangari, A. M. 

2008, AJ, 135, 878 
Meynct, G. & Macdcr, A. 2000, A&A, 361, 101 
Morris, P. W., Eenens, P. R. J., Hanson, M. M., Conti, 

P. S., & Blum, R. D. 1996, ApJ, 470, 597 
Naylor, T. 2009, MNRAS, 399, 432 

Negueruela, L, Gonzalez- Fernandez, C, Marco, A., & 

Clark, J. S. 2011a, A&A, 528, A59-I- 
— . 2011b, A&A, 528, A59-F 

Negueruela, L, Gonzalez-Fernandez, C, Marco, A., Clark, 

J. S., & Martmcz-Nuficz, S. 2010, A&A, 513, A74+ 
Negueruela, L, Marco, A., Herrero, A., & Clark, J. S. 2008, 

A&A, 487, 575 
Nurnberger, D. E. A. & Stankc, T. 2003, A&A, 400, 223 
Palla, F. & Stabler, S. W. 1999, ApJ, 525, 772 
Rathborne, J. M., Brooks, K. J., Burton, M. G., Cohen, 

M., & Bontemps, S. 2004, A&A, 418, 563 
Rauw, G., Crowther, P. A., De Becker, M., Gosset, E., 

Naze, Y., Sana, H., van dcr Hucht, K. A., Vreux, J.-M., 

& WiUiams, P. M. 2005, A&A, 432, 985 
Rieke, G. H. & Lcbofsky, M. J. 1985, ApJ, 288, 618 
Sanders, D. B. & Mirabel, I. F. 1996, ARA&A, 34, 749 
Schaller, G., Schaerer, D., Meynct, G., & Maeder, A. 1992, 

A&AS, 96, 269 
Schnurr, O., Casoh, J., Chene, A.-N., Moffat, A. F. J., & 

St-Louis, N. 2008, MNRAS, 389, L38 
Schnurr, O., Moffat, A. F. J., Villar-Sbaffi, A., St-Louis, 

N., & Morrell, N. L 2009, MNRAS, 395, 823 
Shara, M. M., Moffat, A. F. J., Gcrke, J., Zurck, D., 

Stanonik, K., Doyon, R., Artigau, E., Drissen, L., & 

Villar-Sbaffi, A. 2009, AJ, 138, 402 

Sicss, L., Dufour, E., & Forcstini, M. 2000, A&A, 358, 593 

Smith, N. 2007, AJ, 133, 1034 

Smith, N. & Brooks, K. J. 2007, MNRAS, 379, 1279 
Smith, N., Brooks, K. J., Koribalski, B. S., & Bally, J. 

2006, ApJ, 645, L41 
Smith, N., Povich, M. S., Whitney, B. A., Churchwell, E., 
Bablcr, B. L., Meade, M. R., Bally, J., Gchrz, R. D., Ro- 
bitaillc, T. P., & Stassun, K. G. 2010, MNRAS, 406, 952 
Stead, J. J. & Hoarc, M. G. 2009, MNRAS, 400, 731 
Sung, H., BcsscU, M. S., & Lee, S.-W. 1998, AJ, 115, 734 
Urquhart, J. S., Busfield, A. L., Hoarc, M. G., Lumsdcn, 
S. L., Oudmaijer, R. D., Moore, T. J. T., Gibb, A. G., 
Purcell, C. R., Burton, M. G., & Marechal, L. J. L. 2007, 
A&A, 474, 891 



18 B. Davits et al. 



Urquhart, J. S., Hoare, M. G., Lumsden, S. L., Oudmaijer, 
R. D., Moore, T. J. T., Brook, P. R., Mottram, J. C, 
Davies, B., & Stead, J. J. 2009, A&A, 507, 795 
van der Hucht, K. A. 2001, New A Rev., 45, 135 
Walborn, N. R. & Blades, J. C. 1997, ApJS, 112, 457 
Walborn, N. R., Mai'z-Apellaniz, J., & Barba, R. H. 2002, 

AJ, 124, 1601 
Wallace, L. & Hinkle, K. 1997, ApJS, 111, 445 
Zapata, L. A., Ho, R T. R, Schilke, P., Rodriguez, L. F., 
Menten, K., Palau, A., & Garrod, R. T. 2009, ApJ, 698, 
1422 

Zinnecker, H. & Yorke, H. W. 2007, ARA&A, 45, 481 



